Caging and mosaic lengthscales in plaquette spin models of glasses 
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We consider two systems of Ising spins with plaquette interactions. They are simple models of 
glasses which have dual representations as kinetically constrained systems. These models allow an 
explicit analysis using the mosaic, or entropic droplet, approach of the random first-order transition 
theory of the glass transition. We show that the low temperature states of these systems resemble 
glassy mosaic states, despite the fact that excitations are localized and that there are no static 
singularities. By means of finite size thermodynamics we study a generalised caging effect whereby 
the system is frozen on short lengthscales, but free at larger lengthscales. We find that the freezing 
lengthscales obtained from statics coincide with those relevant to dynamic correlations, as expected 
in the mosaic view. The simple nucleation arguments of the mosaic approach, however, do not 
give the correct relation between freezing lengths and relaxation times, as they do not capture the 
transition states for relaxation. We discuss how these results make a connection between the mosaic 
and the dynamic facilitation views of glass formers. 

PACS numbers: 64.60.Pf, 75.10.Hk 
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I. INTRODUCTION 

The aim of this paper is to study possible connections 
between two perspectives on the glass transition I*! . One 
is the mosaic, or entropic droplet, view that follows from 
the random first-order transition theory of Kirkpatrick, 
Thirumalai and Wolynes (KTW) |1I3: a deeply super- 
cooled liquid is viewed as a patchwork of correlated meso- 
scopic regions which relax by entropically driven nucle- 
ation events. These "mosaics" are characterised by a 
length scale that diverges on approach to a temperature 
Tk, where a Kauzmann entropy crisis occurs. The 
growth of static correlations, under certain assumptions 
for droplet nucleation, leads to a divergence of relaxation 
times of the Vogel-Fulcher kind l] in three dimensions. 

The second, in principle very different, approach is 
based on the idea of dynamic facilitation 0,13: glassi- 
ness is not due to any precipitous thermodynamics, but 
is a consequence of effective constraints in the dynam- 
ics. Here, the central feature is dynamic heterogeneity 
1^0, and the corresponding dynamic scaling, i.e., grow- 
ing times are accompanied by growing dynamic, but not 
necessarily static, correlation lengths [l^- In this view 
there are no finite temperature singularities, and scal- 
ing properties are controlled by zero temperature criti- 
cal points where dynamic lengths and times diverge pJLl| . 
The simplest models that realize this perspective are ki- 
netically constrained lattice models, such as the facili- 
tated spin models of Fredrickson and Andersen (FA) and 
Jackie's East model (6i. il2i..l3i] . 

Here we study lattic e sp in models with plaquette in- 
teractions 0, 0, 0, Il7j from the mosaic or entropic 
droplet perspective. The two models we consider are the 
square plaquette model (SPM) HI . llSll and the trian- 
gular plaquette model (TPM) [l^, [13 • These models 
have exact dual descriptions: one in terms of interacting 
spin variables with standard single-spin Glauber dynam- 
ics; another in terms of free excitations with dynamics 



subject to kinetic constraints. The dynamics of the SPM 
is like that of the facilitated FA model .15 ], th e dynamics 
of the TPM, Ukc that of the East model|17|. 



While plaquette models are realizations of the dynamic 
facilitation scenario, they also allow for a detailed and ex- 
plicit analysis using the mosaic approach. They therefore 
allow a direct analysis of the similarities and differences 
between the facilitation and the mosaic perspectives on 
glasses. In this work, we apply to these models the pro- 
cedure recently suggested by Bouchaud and Biroli (BB) 
[l8l |. by which the mosaic lengthscale is estimated from 
the partition function of finite droplets within a larger 
system. We show that low temperature states do resem- 
ble mosaic states, and that the associated lengths can be 
extracted using the BB procedure. We also show that 
these thermodynamic lengthscales have the same scaling 
as those extracted from bulk many-point static correla- 
tions, and from multi-point dynamical correlation func- 
tions which measure dynamic heterogeneity. In the case 
of the SPM, the static analysis gives two typical caging 
lengths, a feature not anticipated in BB. Moreover, the 
relation between caging lengths and relaxation times in 
these models is not that expected from the mosaic ap- 
proach, as the nucleation assumptions miss the relevant 
transition states and overestimate the free-energy barri- 
ers to droplet melting. 



The paper is organized as follows. In section II we give 
a summary of the arguments of KTW and BB, and make 
general observations about the applicability of this ap- 
proach to the models studied in this work. In section III 
we apply the mosaic procedure to the SPM, and compare 
the lengthscales obtained with other relevant static and 
dynamic lengths of the system. In section IV we repeat 
the analysis for the TPM. Finally, in section V we give 
our conclusions. 
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II. BACKGROUND 
A. Mosaics and their lengthscales 

The basic arguments underlying the BB procedure to 
compute a mosaic lengthscale are the following ^, 18]. 

Imagine a glassy state in which an atom is prevented 
from moving by the fact that its neighbours are fixed 
in position. This picture is consistent only insofar as 
the neighbours of the original atom are fixed by their 
neighbours, and so on. Consider now the original atom 
and its neighbours as a small droplet. One can com- 
pute the probability that the droplet can relax, assuming 
that its boundaries are fixed. If this probability is small 
then rearrangements require correlated motion over dis- 
tances larger than this droplet. One must then increase 
the droplet size and recalculate the probability of rear- 
rangement. In a system with finite-ranged interactions, 
as the droplet is made larger, we eventually cross over to 
a regime where the rearrangement probability is large. 
All larger droplets will be then able to rearrange by se- 
quential moves of this type. 

This crossover defines the mosaic lengthscale 0, Il8| . 
The result is a picture of a glassy state in which droplets 
of small sizes are jammed, or frozen, by their boundaries, 
whereas those with large sizes are unjammed, or melted. 
As the temperature decreases, the typical lengthscale 
separating jammed and unjammed droplets increases. 
(The idea of using finite size scaling to extract a length- 
scale is not at all new, of course. See Refs. 0, for 
applications to glass formers.) 

The above remarks can be made quantitative by way 
of three assumptions 0, E| : 

(i) A droplet of size ^ contains very many metastable 
states whose number scales as e*''^*^'^''^ , where Sc is the 
configurational entropy density which vanishes at the 
Kauzmann temperature, Tr-. 

(ii) The probability of finding a droplet in a given state 
i, given that its boundary conditions are in the state a is 



-/3(?''/.+«''Tfi„) 



(1) 



where 2^ is a partition sum, /; is the free energy density 
of state i and £,^Tia is the free energy cost of matching 
the bulk state i to the boundary state a. The exponent 
9 should satisfy < 6* < d in a short ranged system. 

(iii) Moves involving co-operative motion over a length- 
scale ^ occur on a timescale 



(2) 



where A is a microscopic energy scale, and ip is some ex- 
ponent. If we assume that these moves involve nucleation 
events of one state in a background of an uncorrelated 
state then the exponent V is expected to obey iIj > 6. 

The lengthscale of the mosaic state can be extracted 
by comparing the probability of a droplet state satisfying 
a frozen boundary condition, with the total probability 



of all other droplet states with the same frozen boundary 
[l^. That is, choosing a state i such that Tia — 0, we 
define the crossover length by 



(3) 



The right hand side is a sum over e''''^* terms, so we can 
approximate it by 

g-OMi ^ g-{C?[/3/(T)-3e(T)]+C«/3T(T)}^ ^^-j 

where the sum over states j was replaced with the typ- 
ical number of terms in the sum, *<:(^)^ multiplied by 
the typical weight e~^^^ f(T)+i t(T)]^ have included 
explicitly the temperature dependence of /, Sc and T, 
as they have been averaged over internal configurations 
of the droplet, subject to the boundary condition a. As 
long as the free energy of the well-matched state is typi- 
cal, fi ~ f{T), the mosaic lengthscale is 



(3T 



(5) 



The configurational energy of the droplet is thought to 
vanish at some finite Tr-, which leads to the divergence 
of ^* at that temperature ,2]. 

Combining these thermodynamical arguments, with 
assumption (iii) above we get the typical relaxation time 
for a droplet of size ^* , 



exp(/3Ae 



■4>\ 



(6) 



Droplets of length smaller than will not relax at all in 
this timescale. Droplets of length larger than relax by 
combinations of moves over lengthscales of the order of 
in times of the order of t. Hence, r g ives the typical 
rate of rearrangements in the system j2l| . 



B. Plaquette models 

The BB procedure outlined in the previous subsection 
is extremely general. However, analysing the partition 
function of any model in sufficient detail to evaluate the 
length as defined in Eq. (PJ, is rather difficult. In the 
rest of this work we apply the procedure to two models of 
Ising spins with plaquette interactions in two dimensions. 

We will define the models in the next section, but we 
first make some introductory remarks. The first model 
is the square plaquette model (SPM), in which spins are 
defined on a square lattice, and interactions involve the 
quartets of spins that form the square plaquettes of the 
lattice In the second model, the spins occupy the 

sites of a triangular lattice, and interactions are between 
trip lets of spins on upward pointing triangular plaquettes 
Iff, Tf] . We refer to it as the triangular plaquette model 
(TPM). Both models have dual representations [Tsllr^: 



in a spin representation they describe interacting spins 
with simple dynamical rules (single spin-flips); in a "de- 
fect" representation they describe thermodynamically in- 
dependent excitations with more complicated dynamical 
rules. In this latter description they resemble the ki- 
netically constrained models advocated as effective mod- 
els for the glass transition in the dynamic facilitation 
approach 10]. In this sense, they interpolate between 
models with microscopic degrees of freedom, and more 
phenomenological models of "mobility fields" . 

The SPM and TPM have trivial thermodynamics: the 
excitations in the models are localized and statically non- 
interacting, so the only thermodynamic singularity is at 
zero temperature. In terms of our assumption (i), above, 
this means that we must have Tk = 0. It also has im- 
plications on the assumptions of (ii) and (iii) about the 
surface tension and barrier exponents 9 and ip- Below 
we show that the BB procedure can still be applied to 
extract caging lengthscales from finite-size static proper- 
ties. However, while relaxation times grow with growing 
mosaic lengths, the quantitative dependence of assump- 
tion (iii) does not hold in the SPM and TPM, and these 
relations need to be generalized. 



III. SQUARE PLAQUETTE MODEL 

The SPM is defined by the Hamiltonian [ll[l5l|: 

H = —(1/2) {<yxy<7x,y+lO'x+l,y(^x+l,y + l " 1) , (7) 



where the {fJxy} are Ising spins, a = ±1, and {x^y) indi- 
cates position on a square lattice. The dynamical moves 
are single spin-flips with Glauber rates. If we write 



Pxy 



"^xy^x.y^ l(^x + l,y<^x+l,y+li 



(8) 



then the partition function is simply a sum over the non- 
interacting variables Pxy (up to a non-extensive set of 
constraints on the Pxy from the boundary conditions on 
the spins). The variables Pxy are defined on the pla- 
quettes of the square lattice, which forms a dual square 
lattice, see Fig. ^ Plaquettes with p^y = — 1 cost an 
energy of unity, and are sparse at low temperatures: we 
refer to them as defects. 

The plaquette model has activated dynamics at low 
temperatures. That is, there are very few (if any) possi- 
ble transitions out of a typical state into states with lower 
energy. Thus, low temperature states are very close to 
inherent structures (on quenching to zero temperature, a 
few local relaxations will take place, quickly leading to a 
state with no available transitions). Thus it is an ideal 
model to investigate the mosaic hypothesis: there are 
many metastable states, and these can be easily identi- 
fied. (In contrast, the statistics of the free energy minima 
of an atomistic liquid or a spin glass are much harder to 
probe). 
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FIG. 1: Sketch showing spins (up or down) on vertices of a 
square lattice, and defect variables (circled or blank) on the 
dual lattice formed by the square plaquettes. 



The trivial bulk thermodynamics of the SPM are af- 
fected by freezing the boundary spins of a finite droplet. 
Consider the partition function of a finite system of 
spins (assumed square for convenience), with (4L — 4) 
frozen spins along the boundary. The state of the sys- 
tem is then defined by the configuration of the [L — 2)^ 
remaining (bulk) spins, and its energy is given by (half 
of) the sum of the {L — 1)^ defect variables, Pxy We see 
immediately that there are more configurations of the 
plaquettes than of the bulk spins: a partition sum over 
the bulk spins does not contain all configurations of the 
plaquette variables. Rather, there are 



^plaqs 



^spins 



= {L-lf - {L-2f = 2L-i, (9) 



constraints on the possible arrangements of the defects. 

To understand the origin of these constraints, observe 
that flipping any of the {L ~ 2)^ free spins preserves the 
parity of the number of defects in any row or column of 
the square (dual) lattice. Therefore the 2L—3 constraints 
select the parity of the L—1 rows and L—1 columns of the 
dual lattice. [Fixing the parities of all the rows sets the 
parity of the total number of defects. Thus there are only 
L ~ 2 independent column parities, the final one being 
fixed by the parity of the total number of defects. The 
result is that there are 2L — 3 independent constraints 
on the plaquettes, as required.] These arguments are 
illustrated in Fig.|21 

If we separate the energy of a configuration into a part 
coming from the (L — 3)^ plaquettes in the bulk, and the 
4L — 8 plaquettes along the boundary then we may write 
the partition function for all spins as 



7. — ^-/3(ni+ni, 



(10) 



where i labels the configuration of the bulk spins and a 
that of the boundary. The integer rii is the number of 
bulk plaquettes that are excited, and n^Q, is the number 
of boundary plaquettes that are excited. Making contact 
with section m we identify rii with the free energy fi 
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FIG. 2; A droplet of size L = 7. There are (4L - 4) = 24 
frozen boundary spins (black circles) and (L — 2)^^ = 25 free 
bulk spins (white circles). The energy is determined by the 
state of the 36 defect variables that sit on the dual lattice. 
This energy can be divided into a boundary contribution that 
comes from the sites of the dual lattice marked with +, and a 
bulk contribution from the remaining sites. The partition sum 
is over the 2^^ configurations of the bulk spins; these corre- 
spond to 2^^ of the 2"^^ configurations of the defect variables. 
A given configuration of the defect variables contributes to 
the partition sum if and only if it satisfies the 9 independent 
constraints set by the frozen boundary spins. For example, 
one such constraint is that the parity of the number of de- 
fects in the third row of the dual lattice is the same as the 
frozen spin combination (7i(T2cr3(T4, where the spins (ti_4 are 
identified in the figure. 



and with the surface energy Tia- Note that we have 
replaced the sum in (|3Jl over inherent structures by a sum 
over states, since the typical droplets under consideration 
cannot reduce their energy by any single spin-flip. 

The partition function in can be written in terms 
of partial sums over boundary conditions a. 



(11) 



Now imagine choosing a droplet in an infinite system. In 
doing so we find that its boundary state is a. The droplet 
lengthscale is calculated by considering the contribution 
of that droplet to Za. In other words, we need to estimate 
the expectation value of Zia/Za as a function of droplet 
size. This expectation value is 



Za 



Za 



(12) 



The second form for the expression above emphasises 
that it is an average over the boundary conditions with 
their thermal weights, Za- 

We define ^* as the system size at which the above ex- 
pectation value ceases to be dominated by a single state. 



Since the droplet entropy is 

Sa = -fcs ^(Zia/Za) log{z^a/Za), (14) 
i 

Eq. will be satisfied when (Sa) be of the order of 
ks- In what follows we set /cs = 1. The entropy l(T^ 
measures the typical number of states contributing to the 
partition sum with a fixed boundary. It must not be con- 
fused with the configurational entropy which measures 
the number of states contributing to the overall partition 



A. Finite size thermodynamics 

The partition function for a finite SPM with periodic 
boundaries was calculated by Espriu and Prats ^2^. In 
appendix ^ we generalise their argument to allow for a 
given state of 4(L — 1) fixed boundary spins. The deriva- 
tion is in the spirit of a high temperature series, which 
can be resumed exactly for this model since the bulk 
thermodynamics are those of a trivial free lattice gas. As 
mentioned above, the boundary conditions on the spins 
constrain the parity of the number of defects in each row 
and column of the dual lattice. If the number of rows 
with odd parity is r and the number of columns with 
odd parity is r', then Za depends on the boundary con- 
ditions only through r and r' . Since it is a symmetric 
function of these two integers, it is convenient to define 



m = max(r, r'). 



(15) 



so that Za = Zran- Notc that n is always even since r 
and r' (and therefore m) all have the parity of the total 
number of defects in the droplet. It is also convenient to 
parameterize temperature by: 



(16) 



{Zia/ Z, 



(1/2). 



(13) 



which is related to the average concentration of excita- 
tions by ((1 — pxy)/2) = c/ {1 -\- c) ~ c, the last approxi- 
mate equality being valid at low temperatures. 

The low temperature behaviour of the plaquette model 
obeys scaling relations [2^. Length and timescales di- 
verge at low temperatures as simple powers of c. We 
therefore work in the scaling regime, at leading order in 
c. In terms of droplet sizes, this means that numeri- 
cal values of L may be large. The relevant indicators 
of droplet size will be cL and cL^. That is, we assume 
c ^ 1 ^ L, but make no assumption on the absolute 
sizes of cL or cL^. 

We begin by considering fixed L and very low temper- 
atures: cL^ <^ 1. On picking a droplet from an infinite 
system, it is very likely that it contains no defects. The 
boundaries will then be such that all rows and columns 
contain even numbers of defects. There are 2^^~^ such 
boundary conditions (one for each choice of spins, say, in 
the bottom and leftmost sides of the box). Tracing over 
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FIG. 3: (Top) Droplet entropy density Sqq/L^ versus c = 
for droplets of different sizes L, with frozen boundary 
conditions enforcing even numbers of defects in every row or 
column, m — n = Q. The symbols are finite L entropies 
calculated using the expressions of appendix^ The entropy 
is extensive for cL > 1. On approaching cL ~ 1 from above 
it deviates from the bulk entropy (solid line, L oo), and 
becomes subextensive for cL < 1. (Bottom) Plot of total 
entropy ^oo against cL showing the (non-extensive) scaling 
predicted in H18|l . 
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FIG. 4: Ground states of droplets {L = 6) with different 
boundary conditions. Spins are on vertices on the square 
lattice and are not shown. Excited plaquettes are marked with 
X symbols, (a) Boundary conditions enforce two rows and two 
columns with odd numbers of defects, so that m = 2, n — 0. 
There are two ground states that differ by fiipping the two 
spins identified with black circles, (b) Boundary conditions 
enforce two rows with odd numbers of defects, but all columns 
with even numbers: m = 2, n — 2. There are (L — 1) = 5 
distinct ground ground states of which two are shown: they 
differ only at the four marked spins 



the bulk spins with each configuration of the boundary 
contributes Za = Zqq = •^m=o,ri=o to Z: 

Z ^ 2'-'^-'[Zoo + 0{cL^)]. (17) 

Now, if m = n = then the boundary of the droplet 
specifies its unic^ue ground state. Further, simple count- 
ing arguments show that excited states contain even 
numbers of defects u > 4. Their Boltzmann weights are 
c" and the degeneracy of the state with u defects is pro- 
portional to L". This means that: 

Zoo = l + OicL). (18) 

This partition function is also relevant to periodic bound- 
ary conditions, and is the one calculated in In Fig. 
iniwe show the droplet entropy 6*00 = dTiTln Zqq), as a 
function of c for different droplet sizes L. The top panel 
of Fig. O shows that the droplet has an "entropy crisis" 
at cL ^ 1 . In the bottom panel of Fig. O we see that 
the droplet entropy scales with cL when cL <^ 1, and is 
subextensive in this regime. 

From figure |31 it is clear that as the temperature is 
increased, droplets with m = n = remain frozen un- 



til cL ^ I. However, Zqq dominates the droplet parti- 
tion function only for cL^ 1. This means that typical 
droplets are frozen only for cL^ ^ 1. Between cL ^ 1 
and cL^ ~ 1 we must consider Z,„„ with to, n 7^ 0. 

The behaviour of droplets changes qualitatively in the 
case of ground states containing two or more defects. As 
is clear from Fig.^ if m > 2 then there is more than one 
state of the bulk spins that satisfies the boundary condi- 
tion perfectly. Specifically, if m = 2 then there are two 
ground states, each of which contains two defects, and 
contributes to the partition function. More generally, 
we find that 

Z„o-m!c'"[l + 0(cL)]. (19) 

It follows that the zero temperature droplet entropy is 
now Smo{T = 0) = log to!, and so for these boundary 
states, we expect typical values of 

z^c./Zo, ~ (1/m!) + 0{cL). (20) 

While the droplets do have a feature in the entropy 
at cL 1, the degenerate ground states all contribute 
equally to the partition sum at low temperature. Thus, 
as soon as droplets with m > 2 dominate Z then the 
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FIG. 5: Entropy as a function of cL for L = 100. The full 
line is the bulk entropy of the SPM. The open symbols are the 
entropies SmO for frozen boundaries with different m, as ob- 
tained from the finite size partition functions of appendix 1X1 
The + symbols give the entropy in the case where boundary 
conditions are generated randomly. 
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they are thermodynamically melted according to condi- 
tion H13|l . Fig. [SI shows that this happens as soon as 
cLF' ~ 1, as expected. The melting lengthscale is there- 
fore, 



-1/2 



(21) 



It is clear that there are two crossovers in the droplet 
thermodynamics as c is increased. For cL^ ^ 1 nearly 
all droplets are frozen in true ground states of the sys- 
tem. For cL ^ 1 the system behaves as the bulk. In the 
intermediate regime cL ^ 1 <C cL^, each droplet is al- 
most certainly in a state that minimises the energy with 
respect to the boundaries; however, there are many of 
these states for a fixed configuration of the boundary. In 
this regime m is about cL^ . The low temperature droplet 
entropy is then 



St= 



cL^ log cL^ 



(22) 



which is large compared to unity, but small compared to 
its bulk value cL^ log 1 /c (since cL <^ 1 implies cL^ <C 
1/c). We show this explicitly in Fig. [S] where we average 
the exact results for the entropy at given m and n over 
a distribution of mn that is sampled by Monte Carlo. 

Thus far we have analyzed only states with n — 0. 
That is, the boundary condition requires m rows to have 
odd numbers of defects in them, and it also requires m 
columns to have the same property. The probability that 
?z = is the probability that no two defects occupy the 
same row or columns as each other, 

P{n = 0) = 1 - 0{c^L^). (23) 

This probability is close to 1 as long as cL^ ^ VL- 

The melting length H21(l can also be accessed through 
static overlap functions, which will be important in our 
discussion of dynamics. Consider the self-overlap for a 
given boundary a, 



Coa — 



1 



{L-2y 



(24) 
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FIG. 6: Plots of the overlap function {axy)a, averaged over 
boundary conditions a at various temperatures. (Top) Aver- 
age overlap Co as a function of temperature. (Middle) Col- 
lapse of this data as a function of cL^ (the dashed line is a 
guide to the eye) . (Bottom) Spatial dependence along the di- 
agonal of the square droplet, showing strong correlations near 
the boundary. 



where we have explicitly summed over the spins of a 
droplet with given boundary condition, since for a given 
realisation of the a, single spin correlations are finite and 
depend on position. The brackets (. . .)a denote a ther- 
mal average over configurations of the bulk spins, with 
the boundary spins fixed. Coa should give the long time 
limit of the spin-spin autocorrelation for a droplet with 
given boundary conditions. Averaging over the boundary 
spins we get: 



Co 



(25) 



xy 



which is analogous to a disorder average. 

In Fig.|H|we show these overlaps. The correlator {(Jxy)a 
can be calculated in a high temperature series: see ap- 
pendix ^ The top and middle panels of Fig. show 
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that Co is a function of cL^ , which is close to 1 for small 
cL^, and becomes small when cL^ gets large. This scal- 
ing also indicates that Co will be small in the regime 
cL <C 1 ^ cL^. The bottom panel of Fig. |H| shows the 
spatial dependence of {iyxy)a- The spins near the bound- 
aries are strongly fixed by the boundary condition. This 
correlation decays away from the boundary with a length- 
scale that scales as ^4 ~ c~^/^. These observations for 
the overlaps are consistent with thermodynamic melting 
of the droplet at cL^ ^ 1 found above. 

B. SPM and mosaic lengths 

We have established that the melting lengthscale ^* 
scales as c~^/^ in the SPM, but that there is also a second 
important thermodynamic crossover at L ^ c^^ > 
We now compare the results so far with the expectations 
of BB. The appearance of two lengthscales was not an- 
ticipated in [l^. We have shown that the crossover at 
cL^ ~ 1 comes from a small number of states that are 
have similar weights in the droplet partition function. 
We find that distinguishing these few specific states from 
typical states with bulk energy E is necessary when dis- 
cussing the droplet partition function Z^. 

We use to define the mosaic lengthscale The 
assumption of BB is that the lengthscale defined by Q 
will scale in a similar way. Eq. Q treats all states in- 
dividually, whereas Q assumes that we can extrapolate 
the properties of the partition sum from typical states at 
the relevant temperature. We write the mosaic length- 
scale extracted from Q as C*,typ: we find that and 
^*,typ do not coincide in the SPM. 

More precisely, the average bulk energy of a droplet at 
a given temperature is c{L — 2)^, and the configurational 
entropy at that temperature is c{L — 2)^ ln(l/c) -I- {2L — 
l)ln2. If we replace the bulk of the droplet by a new 
bulk state with the same energy, then the typical bound- 
ary energy is 2(L — 2) (half of the boundary plaquettes 
will be excited, since there are no two-point spin correla- 
tions either in the bulk or along the boundary). So from 
Eq. (gj, we have = 2, 6* = 1 and 

On the other hand, we showed in the previous section 
that the Eq. © leads to 

The origin of this discrepancy is simply the degenerate 
ground states of the droplet, whose boundary energy does 
not take the typical value T(r) ~ L, but rather have 
Tia = 0. Thus they must be separated from the sum 
before it is approximated: 

+ J2 e-^hL^-P^.^ (26) 



where the sum over j G gs denotes a sum over states 
with — 0. The final term in H26|) can then be approx- 
imated to give a form resembling Eq. (0J, but that term is 
irrelevant to droplet melting in the SPM since the second 
term dominates the partition sum for cL <C 1 ^ cl? . 

We interpret this result as evidence that approxima- 
tions such as those leading to Eq. Q are rather danger- 
ous. In general, systems may possess significant numbers 
of states in the partition sum with smaller than aver- 
age boundary energy, and these may be sufficient to de- 
stroy spin correlations in the infinite time limit (Co — > 
when L <C C*,typ)i even if the droplet is exploring only a 
small fraction of its metastable states. Further, we will 
show in the next section that if the droplet does melt in 
stages then choosing whether to identify the bulk relax- 
ation time as t(^,) or T(^,_typ) requires investigation of 
the specific system under study. In the SPM, we will find 
that while complete destruction of of spin correlations at 
large times does take place for L ^ C*,typj the timescale 
for this relaxation is much longer than the bulk relax- 
ation time. However, this effect cannot be inferred from 
thermodynamic arguments: it is a kinetic effect. 

To end this discussion, we note that while the condi- 
tion cl? ~ 1 leads to rather small droplets, the configura- 
tional entropy of a droplet of this size is cl? log(l/c) ^ 1. 
BB comment that the vanishing of the configurational en- 
tropy leads to a trivial lower bound on the mosaic length- 
scale which in this case is: Q > [clog(l/c)]^^. Thus we 
see that despite the condition that frozen droplets con- 
tain at most one defect, the freezing lengthscale is still 
much larger than its trivial lower bound. 



C. Relaxational dynamics 

In this subsection we consider the relaxation of the 
droplets studied thermodynamically above. The relevant 
correlation function is 

^^«W=(^^ E MOMt))o^ (27) 
a; ye bulk 

where the average is over realisations of the Markov chain 
that describes the time evolution of steady state of the 
system (with fixed boundary state a) . We also define the 
average of Ca{t) over boundary conditions 

C{t)^Y^Mt) (28) 

a 

The long time limit of Ca(t) is given by the overlap Coq. 

^From Monte Carlo simulations, it is clear that the 
relaxation of the system slows down as the system size is 
decreased towards 1/c: see Fig. [7| From our discussions 
of the static single point correlator [Eq. (|24|l and Fig. |H] 
we expect the long time limit of this function to approach 
zero even for cL <C 1 <C cL^, although this is hard to 
show in simulations since simulations with large numbers 
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FIG. 7: Simulations of the SPM at c = 0.01 for various sys- 
tem sizes and frozen boundary conditions. The relaxation 
time increases significantly as cL decreases through unity, but 
the long time limit of C{t) would still appear to vanish since 
cL^ ^ 1. We also show relaxation for L = 1000 and periodic 
boundary conditions. The relaxation can be well-fitted by an 
exponential for large L. At smaller L the decay is slower than 
exponential. 
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FIG. 8: Simulations of the SPM at c = 0.01 for various system 
sizes and frozen boundary conditions. The long time limit of 
C{t) increases from zero as cL^ becomes of order unity. 



of defects at small temperatures are very computationally 
expensive. 

Moving towards the regime with cL^ '--^ 1, we see that 
the relaxation of the system is very slow, and that there 
are large correlations that persist to the infinite time 
limit: see Fig. |H1 We can make a link between static 
and dynamic quantities by evaluation of the long time 
limit of these functions using the high temperature series 
(cf. Fig. 0). 

Thus we have shown how the two lengthscales asso- 
ciated with the glassy mosaic state affect the dynamics 
of the system. As the temperature is reduced, the re- 
laxation slows for cL ~ 1, before stopping completely at 
cL^ ~ 1. The sep aration of these lengthscales was not 
anticipated in |l8l |: we have more to say on this below. 

Before ending our discussion of dynamics, we make 
two further points. Firstly, we note that it is clear from 
Ref. ^3 that the bulk relaxation is caused by correlated 
motion of defects over distances ^dyn ~ c~^. This coin- 
cides with the system size above which the system relaxes 
as the bulk. Secondly, it is also known [23, [23| that cor- 



related motion for defects over a distance / takes place 
on a timescale t ~ le'^^ so we see that the analogue for 
equation lO is 

Tr^e^^+'°si (29) 

The exponent depends logarithmically on ^ with a pref- 
actor that does not depend on temperature: apart from 
the constant activation energy, the free energy barriers 
to motion are purely entropic in origin. The entropic na- 
ture of the barriers explains why relaxation of the SPM 
model diverges only in an Arrhenius fashion (r ~ e'^'') 
despite the co-operativity of its dynamics. 



D. SPM summary 

It is clear from the results of this section that, even 
in a simple model like the SPM, the situation regarding 
caging and mosaic lengthscales is already more compli- 
cated than what the appealingly simple and general ar- 
guments of would predict. As we increase the size of 
droplets with frozen boundaries, there are two crossovers. 
Static spin-spin correlations vanish at at c ~ but 
neither the droplet entropy, nor its dynamical correla- 
tion functions reach their bulk forms until the much 
higher temperature c ^ L~^. We note that these two 
lengthscales were also identified in the discussion of static 
and dynamic lengthscales in Ref. [2^. To be precise, 
the static four-spin correlator (o'ooO'orO'roO'rr) decays on 
a lengthscale ^1'*^* ^ c~^^'^, but the dynamic four point 
correlator decays on a length that scales as ^2,2 ~ c~^. 
There is also a static lengthscale that mirrors the scaling 
of ^2,2 which is derived from fluctuations of two point 
correlations [23l |. 

Despite these complications, a general picture like that 
of 18] does seem to be applicable to the SPM, in the 
sense that droplet melting does occur, even if it hap- 
pens in stages. As the temperature decreases towards 
the Kauzmann singularity at T = 0, the mosaic length- 
scale increases, and the relaxation time increases accord- 
ing to H29|l . This increase follows a simple Arrhenius law 
because only the entropic barriers depend on This 
suggests that one should replace |2Jl by the more general 
parametrization 

r ^ fsiOe^^-^^^ (30) 

where /s(^) describes the entropic part of the barrier and 
fuiO the energetic part. At small temperatures then we 
expect the energetic part to dominate, except in the case 
where fu{0 independent of 1^, as is the case for the 
SPM. 

It is at this point at which our discussion of dynam- 
ics diverges from that of KTW and BB: since the (non- 
perturbative) dynamical moves are not nucleation events 
in this model, then the free energy barrier is not set by a 
transition state with a critical droplet of the new phase. 
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FIG. 9: Sketch showing relation between square and rhom- 
boid droplets {L = 6). The spins are on the intersections 
of grids; an is marked by a filled circle. The x symbol 
marks the plaquette interaction §24 ~ cr240"25cr34. Rotations 
of 120° leave both Hamiltonian and triangular lattices invari- 
ant. However, note that the boundaries of the droplets are 
not all equivalent, since the rhombus shape is not invariant 
under this rotation. 



Rather, the relaxation proceeds by less expensive rear- 
rangements that occur over a lengthscale ^*,typ, as was 
discussed in detail in Ref . 1231 . 




FIG. 10: Illustration of Eq. ^ with k ^ 6. The product of 
the five spins marked by filled circles is given by the product of 
the fifteen plaquette variables marked with crosses. Both the 
marked defects and the spins marked with circles (filled and 
empty) form parts of Sierpinski triangles: they are marked if 
and only if the relevant entry of Pascal's triangle is odd. 



IV. TRIANGULAR PLAQUETTE MODEL 

In this section we compare the results we found above 
for the SPM, with a similar model of Ising spins: the tri- 
angular plaquette model, or TPM 0, llTl . This model 
has a super- Arrhenius divergence of the relaxation time 
at low temperatures, and is therefore a fragile model, 
as compared to the SPM, which is strong. Moreover, dy- 
namical correlations do not have the strong anisotropy of 
the SPM "^^l . The calculation of the partition function, 
however, is considerably more involved, so our results 
below are less detailed than those for the SPM. 

The TPM is defined by the Hamiltonian [HIS: 

H = -{l/2)''^axy(Jx+i,yO'x.,v+i- (31) 

xy 

where as before {(Jxy} are Ising spins on a square lattice 
and the dynamics are is spin-flips with Glauber rates. If 
we define the plaquette variables 

Qxy ^xy^x+l,y^x,y + l 

(32) 

we again find that the partition sum over the spins re- 
duces to a product over q^y, that are independent up 
to a non-extensive set of constraints from the boundary 
conditions p^lT^. 

We note that defining this model in this way leads to 
a Hamiltonian that does not possess the square symme- 
try group of the lattice: it is more intuitive to deform 
the square lattice into a triangular one. The interactions 
are then around upward pointing triangular plaquettes 
(see Fig. I^l. The symmetries of the Hamiltonian are 
then translations by lattice vectors, and 120° rotations 
(60° rotations leave the lattice invariant but flip down- 
ward pointing triangles to upward pointing ones, so they 
change the Hamiltonian). 

The analogy of the square droplet of Fig.His a rhombus 
shape in the triangular representation. The arguments 



for counting bulk and boundary states are the same as 
those of the plaquette model. However, identifying the 
22L-1 degeneracies of each state and the nature of the 
(2L — 3) constraints is considerably more involved. 

A. Static spin-spin correlations and constraints 

The key to understanding static properties of the tri- 
angular plaquette model is to realise that 

{q^^,n1^2,l2 ■ ■ ■ <l^r.,u) = [tanh(/3/2)] " . (33) 

All multi-spin correlators that can be written in this form 
are finite. Further, all other multi-spin correlators are 
zero . While one and two point static spin correlations 
vanish, 

{<^xy) = 0, {(Txy(yx'y') = S^^x'Sy^y', (34) 

specific three point correlators are non-zero, 

{axyax+i,y<Jx,y+i) = (qxy) = tanh(/3/2). (35) 
A useful general relation is that 

fe k k'-l 

'^^^yYl'^x+r,y-k= n Hlx+r^y-k'^ (36) 
r=0 k' = l r=0 

where the Bnr = nl/[rl(n — r)!] are binomial coefficients. 
Since the g's and the cr's are Ising variables, only the 
parities of the binomial coefficients are relevant. In terms 
of correlation functions, this relation implies that there is 
a finite multi-spin correlator connecting axy and various 
spins in row a: — fc of the lattice. 

This correlator is related to the first k rows of Pascal's 
triangle (PT). The correlator connects the spins corre- 
sponding to entries of the kth row of the PT that have 
odd parity. Its value is given by tanh(/3/2)^' where p is 
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FIG. 11: Relaxation of finite systems at inverse temperature 
/3 = 4 (c = 0.018). Notice that the long-time limit of C(t) is 
changing, but that the relaxation time is approximately con- 
stant. Compare to Fig.|H] where the relaxation time increases 
sharply as C{t oo) increases. 



the number of entries in the first k— 1 rows of the PT that 
have odd parity. See Fig. ITUland the associated caption. 
Now, the typical number of odd entries in a Pascal's tri- 
angle of linear size I is Z''* with df — log2(3). It follows 
that typical correlations over a lengthscale r scale as 



tanh(^/2)('-''') 



exp 



(37) 



from which we infer that static correlations decay on a 
lengthscale 



6- 



-i/df 



(38) 



In the case of finite droplets there must be (2L — 1) sym- 
metry operations that leave the energy of a droplet in- 
variant, but require the boundary spins to change (recall 
Fig- 121) • Furthermore, fixing the boundary spins imposes 
(2L — 3) constraints on the allowed configurations of the 
plaquette variables. The form of these symmetries and 
constraints are given explicitly in appendix IbI 

The SPM was simple to analyse since the constraints 
are in a form where each plaquette appears in exactly 
two constraints. In that case we can make a geometri- 
cal interpretation of those constraints. In the TPM the 
consequences of imposing the constraints is less clear. 
We therefore proceed directly to a discussion of dynamic 
properties, in an attempt to access the static properties 
indirectly. 



B. Dynamic spin correlations 

Relaxation in the TPM is frustrated by a hierarchy 
of energy barriers (l^ Il7| : co-operative motion over a 
distance ^ requires the crossing of an energy barrier of a 
height log2 ^. This behaviour is reminiscent of the East 
model p^ . 

In order to understand the scaling of the mosaic length, 
we investigate the scaling of C(t) in finite droplets: see 
figure 111! A lower bound on the freezing lengthscale is 



FIG. 12: Relaxation at different temperatures showing bulk 
results (circles) and two different values of cL"^^ . The actual 
system sizes used are (left) L = 13, 26; (middle) L = 16, 32; 
(right) L — 26, 44. The dashed lines are guides to the eye 
showing that the long time limit scales as a function of cL'^' 
only (there are weak deviations from scaling which we at- 
tribute to subleading corrections in c and 1/L). 



set by the elementary condition cL^ ~ 1, since all states 
with single defects have different boundary conditions. 
This would imply that the long time limit of the auto- 
correlation, Co should be a scaling function of cL^, as in 
the SPM, Fig. El 

However, from the discussions of static correlations 
above, it is clear that the physics of the triangular pla- 
quette model is intrinsically related to the parities of bi- 
nomial co-efficients, and hence to Sierpinski's triangle. 
As mentioned above, this structure has a fractal dimen- 
sion of di = log2(3) — 1.585. In Fig. ^1 we show relax- 
ation data at various temperatures with constant values 
of cL'^f . We see that Co does seem to depend only on this 
combination. The conclusion is that the mosaic length- 
scale scales as 



-i/dt 



(39) 



This relation indicates that typical droplets melt when 



they contain of the order of 



defects; be- 



low this number there no co-operative mechanisms allow- 
ing the droplet to rearrange. This should be compared 
with the SPM, in which droplets with L ^ contain 
of the order of one defect, but those at the dynamical 
crossover {cL ^ 1) contain of the order of defects. 

An alternative measure of the co-operativity of the dy- 
namics is given by the bulk "four-point" correlation func- 
tion (see e.g. p^) 



lit)<JOQit)<JxyiO)<JOQ{t)) . 



(40) 



This correlation function measures spatial correlations in 
the single site autocorrelation function. It is convenient 
here to work with the connected part of this function, 
normalised according to 



G4{xy;t) 



1 - ('Too(t)^oo(0))2 



(41) 



In the triangular plaquette model, the function G4 has a 
complicated spatio-temporal structure, again related to 
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FIG. 13: Circular average of the normalised four point dy- 
namic correlator. The system size is L = 128 with periodic 
boundaries, which is large enough to avoid finite size effects. 
(Left panel) Three different temperatures, with times chosen 
so that C{t) — 0.5. (Right panel) Collapse of the data of the 
left panel on rescaling of length as rc^^"^^ . 



FIG. 14: Data showing relaxation time for large systems as 
a function of inverse temperature, p. The relaxation time is 
defined by C(r) = 0.5. The error bars are smaller than the 
symbols shown. The dashed line is a fit to the form Inr = 
ao + aif3 + qf3^ / In 3. Good fits are also possible with different 
coefficients for the quadratic term. 



Sierpinski's triangle: we postpone a detailed discussion 
of this structure to a later work. However, the circular 
average of this function, which we denote by G'4r(r, 
decays smoothly with distance (at a fixed time). This is 
shown in Fig. ^1 from which it is clear that the length- 
scale over which dynamical correlations occur also scales 
as 

We have used the dynamics of finite droplets to identify 
a single mosaic lengthscale in the triangular plaq uette 
model, that scales as c~^^'^' . Further, we know [13 that 
energy barriers to motion over a lengthscale I require an 
activation energy of the order of log2 I. At very low tem- 
peratures, we therefore expect the bulk relaxation time 
at very low temperatures to scale as 

T ~ exp iP log2 ~ exp (P^ / In 3) . (43) 

Note that Ref. predicted exp(/3^/ln4) based on the 
assumption that ^ c^^^. Fig. 1141 shows a fit of the 
relaxation time of C (t) with the function In r = oq + 
ai/3 -|- In 3. The constant and linear terms come from 
entropic and dynamical effects which become irrelevant 
at very low temperatures where the energy barriers on ^ 
dominate, see Eq. (|30ll . The fit in the figure is consistent 
with H43|l . although the range of f3 that we consider is 
too small to confirm it beyond doubt. 



C. TPM summary 

The thermodynamics of finite droplets in the triangu- 
lar plaquette model are not sufficiently simple to allow an 
analysis of the droplet entropy and hence determination 
of the mosaic lengthscale according to the criterion 
However, it appears that the long time limit of the single 
spin autocorrelation function approaches unity as cL^°^^ ^ 



gets small. If Co approaches unity then it seems reason- 
able to assume that droplets are effectively trapped in a 
single state that dominates the partition sum, and that 
this should be an equivalent criterion to 

We find that the mosaic lengthscale extracted in this 
way {£,*), scales in the same way as the static length- 
scale for multi-spin correlations (^stat), and as the dy- 
namical lengthscale extracted from the four point corre- 
lator (^dyn)- This scaling is evidence for the presence of 
a single lengthscale in this model. This is in contrast to 
the TPM, and more in line with 18]. 

On the relation of dynamics to statics, the analogue of 
^ seems to be 

r-fs{Oe^'°^^^ (44) 

The increase of the energy barriers with f is again slower 
than the power law predicted by [l8| , but the assumption 
of an increasing function is robust. The absence of a 
power law dependence again indicates that the system 
relaxes by a more efficient mechanism than nucleation of 
a new droplet in the background of an uncorrelated state. 

V. CONCLUSIONS 

We have successfully applied the procedure of Ref. ^3 
to two simple models of interacting Ising spins. We find 
that the concept of "caging" whereby the states of in- 
dividual spins are strongly constrained by frozen spins 
surrounding them applies well in these models. We were 
able to study the breakdown of this caging effect by con- 
sidering droplets of increasing size. 

In the square plaquette model (SPM), the caging 
breaks down in two stages as the droplet size increases. 
For L ^ c^^/^ the droplet is frozen. For c^^/^ <C 
L <^ cr^ the droplet will eventually relax, but this re- 
laxation is slower than that of the bulk system. Finally, 
for L ^ the relaxation time coincides with that of 
the bulk. The correlation lengths controlling these two 
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crossovers have the same scahngs as the two separate 
lengths identified in this model from four point static and 
dynamic correlations [2^. The two stage droplet melt- 
ing in this model is more complicated than the simple 
scenario of BB. Interestingly, the two crossover tempera- 
tures of the SPM, as well as the structure of the droplet 
states in the intermediate regime, bear some resemblance 
with the situation in mean-field systems like the p-spin 
model (see e.g. and references therein). 

In the triangular plaquette model (TPM) we find only 
a single crossover in the system at L ~ This 
is more in line with the situation anticipated in BB. The 
same lengthscale is again present in static multi-spin cor- 
relation functions, and in four-point dynamical correla- 
tors. 

We draw two conclusions from these results. Firstly, 
we are able to verify that the procedure of BB is a suit- 
able way to identify a mosaic lengthscale. In the general 
case it may be difficult to test the criterion (|3Jl directly, 
although it is possible in the SPM. In less favourable 
situations, it is still possible to use dynamical studies 
to extract the mosaic lengthscale, as we showed for the 
TPM. 

Secondly, we have shown that models without a finite 
temperature phase transition can still be described from 
the mosaic perspective of KTW and BB. This indicates 
that the concept of a glassy mosaic state is more general 
than the specific picture of KTW. While the low temper- 
ature states of the plaquette models are constructed by 
tiling the plane with different metastable configurations 
of the spin system, there are no well-defined droplets or 
domain walls in these states. Rather, there are point- 
like defects that allow the system to interpolate smoothly 
between different metastable configurations. The crucial 
difference with the KTW approach is that, while dynam- 
ics does proceed by non-perturbative events that require 
co-operation over a lengthscale that diverges at low tem- 
peratures, these events are not related to nucleation, but 
to more efficient rearrangements of the point-like defects. 

In summary, typical states in the plaquette models are 
glassy mosaics whose dynamics should be understood in 
terms of point defects with co-operative dynamical rules. 
These defects are not locally conserved, and are not re- 
lated to any concept of free volume. Rather, they rep- 
resent regions in which the free energy barriers to mo- 
tion are smaller than average. In this approach, the 
"activated processes" that dominate glassy a-relaxation 
involve co-operative rearrangements of defects over a 
lengthscale ^. The time scales for these processes increase 
with ^, with a general relation similar to H3U|I . 

The purpose of this work has been to investigate rather 
general arguments of KTW and BB explicitly in specific, 
and rather simple, finite dimensional models. We find 
good evidence for a generalised caging effect, but our view 
of static and dynamic properties of the glassy mosaic is 
qualitatively different from that of KTW. The possible 
links between spatial variation in the mosaic length and 
dynamical heterogeneity remain as an area for further 



study. 
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APPENDIX A: HIGH TEMPERATURE SERIES 
IN THE SPM 

We follow [23 in making a high temperature expansion 
for the square plaquette model. The partition function 
for a droplet of size L with fixed boundary spins is 

Z„ = e-(^-l)'^/2 ^ j-J g-/3p.„/2^ (Al) 

bulk a^y l<xv<L-l 

with 

Pxy — '^xy^x-\-l.y^x,y-\-l^x+l,y+l^ 

(A2) 



as above. Since the Pxy are Ising variables then this can 
be written 




where z tanh(/?/2) = (1 - c)/(l + c). 

Expanding the product results in a power series in z: 
there are 2*^-^"^^ terms corresponding to states of a dual 
system where the spins are the Pxy Each state is equal to 
a power of z multiplied by a combination of the <Txy] the 
term vanishes on summation over the cr^y if it contains 
any spin to an odd power. Otherwise, the coefficient of 
the power of z depends only on boundary spins: in that 
case, the sum over bulk spins just contributes a trivial 
factor 2(^-2)^^ 

Terms containing bulk spins to even powers are formed 
by flipping whole rows and columns of the dual Pxy vari- 
ables. The result is that 

Z„ = 2-2i+3(i + c)(-^-i)'(l/2) X 

^ .^{H+V){L-l)-2HV ^ 

W {ciyCri^y+icrLya-L.y+i)^^ x 

l<y<L-l 

io-xlcrx+l,lcrxLCrx+l,Ly" , (A4) 

l<x<L-l 

where H = "^yhy and V — X]a:^a;i ~ 1 means 
the Pxy in row y have been flipped. We note that the 
boundary spins enter only through combinations such as 
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o'ij/Cri.j/+iO'Lj,crL.j,+i, which gives the parity of the num- 
ber of defects in the yth row of the dual lattice. Fur- 
ther, the terms in tlie sum depend on only four num- 
bers: H; V; the number of rows for which hy = 1 and 
(^iy<^i,y+iO'Ly<yL,y+i = "1; and the number of columns 
for which Vx = I and (Txi(Jx+i,ictxL(Tx+i,l = -1- Let r be 
the number of rows for which aiyai^y+iaLyCrL,y+i = —1 
and r' be the equivalent number of columns. Then we 
can write 

Z„ = 2-2^+2(1 + c)(^-i)' 

r r' L — l — rL — l — r' 

xEE E E <i.>" 



h=Ov=0 h'=0 v'=0 



with 



^h,v,h',v' ~ ^rJiBr' ,vBL-l-r,h'BL-l-r',v' X 

^_l)h+v^(H+V)(L-l)-2HV^ (A6) 

where now H = h + h' and V = v + v' since h and 
h' are the number of rows with hi = 1 and odd and 
even defect parities respectively, and v and v' are similar 
numbers of columns. We denote the binomial coefficients 
by -B„,r = n!/r!(n - r)! 

We may now sum over v and v', yielding 

= 2-2^+2(1 + c)(^-^)' 

r L-l-r 

xE E Br,hBL-l-r,h'Ml% (A7) 
h=0 h'=0 



with 



(A8) 



A similar procedure leads to the expression for the ex- 
pectation value of any spin, given a boundary condition. 
We state only the result which depends on the coordi- 
nates of the spin 2 < xy < L — 1, as, well as the total 
numbers of rows and columns with odd defect parity : r 
and r' . It also depends on s, the mmibcr of such rows 
with y co-ordinate less than j and s': the number of such 
columns with x co-ordinate less than i. (Clearly s < r 
and s' < r'). 

The result is that 



where 



y—l—s r — s L — y + s — r 



(A9) 



^" = E E E E <l,H",H"'^ (Aio) 



with 



h=0 h'=0 h"=0 h"'=0 



^^^h,h',h",h"' — 



^__l^h+h" ^^H ^ ^L-l-H^L-x+s'-r' ^ 

I^^H-2(h+h') ^ ^L-l-H+2{h+h')y-l-s' ^ 
^^H-2(h+h') _ ^L-l-H+2(h+h')y' ^ 

Bs' .hBy-i^s' M' Br' -s' ,h" B^-y+s' -r' M'" (All) 



These results are still a little cumbersome, but they 
are exact, and are can be evaluated computationally. The 

number of terms in the sum may be as large as 0(L*), but 
this is polynomial rather than exponential in the system 
size. 



APPENDIX B: TPM SYMMETRIES AND 
CONSTRAINTS 

Here we make some comments about the thermody- 
namics of finite {L x L) droplets in the triangular pla- 
quette model. We first identify the 2^^~^ transforma- 
tions that leave the droplet energy unchanged. Recall 
that 



'^x+ljy — Q.xy^xy^x,y-\-\' 



(Bl) 



It follows that if we specify the boundary spins cjiy and 
<^xL (where x,y = 1 . . . L) then the rest of the spins are 
determined by the set {qxy} (with x,y = l..._L — 1). 
Thus, the required degenerate states for a given choice 
of the {qxy} can be generated by cycling through the 
possible states of these two boundaries. 

To identify the constraints on the plaquette variables 
arising from the boundary conditions, observe that for a 
groimd state {qxy = 1 V xy) with specified spins along two 
boundaries ({ciy, CTa;/,}) then the state of the remaining 
boundary spins are 



^Si = Ui^iyY'-''^-^ (B2) 

< = n(^i.')^— 'n(^xL)^— ,(B3) 



x=l 



where a^"'" is equal to unity unless a = —1 and J3„,r is 
odd, in which case it takes the value —1. 

For a general configuration of the g''s and the boundary 
spins oiy^GxL, we have 



(B4) 



x—l x—y 

<yxl = a'x.Yl Hiqx'yf-^'-^-y-^, 
j/=l x' = l 

L-y-1 L-y'-l 

^Ly = n n iiv,v+v'f'^—''''- (B5) 



=0 x=l 



Multiplying both sides of each constraint by the appro- 
priate cr" leaves an equality between a product of bound- 
ary spins and a product of plaquette variable. These are 
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the 2L — 3 constraints that determine which configura- 
tions of the excited plaquettes are present in the partition 
sum. 

While these expressions appear comphcated, for a 
given droplet boundary condition they are simply con- 
straints of the form 



±1, 



(B6) 



In the SPM the constraints were of the form 



L-l 

YIpccv =±1. 



(B8) 



where the constrained plaquettes form parts of Sierpinski 
triangles. For example, for L = A, the constraint coming 
from the boundary spin (742 is 



9l2922932<Z23 — Cri2Cri3Cr240'42 ■ 



(B7) 



The simplicity of (|B8|I in the SPM, as compared to (|B5|I . 
is the reason for the relative intractability of the partition 
function of the TPM. 
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